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In a multi-machine environment, the inter-area low-frequency oscillations 
induced due to small perturbation(s) has a significant adverse effect on the 
maximum limit of power transfer capacity of power system. Conventionally, 
to address this issue, power systems were equipped with lead-lag power 
system stabilizers (CPSS) for damping oscillations of low-frequency. In 
recent years the research was directed towards optimal control theory to 
design an optimal linear-quadratic-regultor (LQR) for stabilizing power 
system against the small perturbation(s). The optimal control theory provides 
a systematic way to design an optimal LQR with sufficient stability margins. 
Hence, LQR provides an improved level of performance than CPSS over 
broad-range of operating conditions. The process of designing of optimal 
LQR involves optimization of associated state (Q) and control (R) weights. 
This paper presents an analytical approach (AA) to design an optimal LQR 
by deriving algebraic equations for evaluating optimal elements for weight 
matrix ‘Q’. The performance of the proposed LQR is studied on an JEEE test 
system comprising 4-generators and 10-busbars. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Sreenivas Uravakonda 


Department of Electrical and Electronics Engineering, Faculty of Engineering 

Srinivasa Ramanujan Institute of Technology, Rotarypuram, Jawaharlal Nehru Technological University 
Anantapur, Ananthapuramu, Andhra Pradesh, India 

Email: uravakonda.sreenivas @gmail.com 


1. INTRODUCTION 


To enhance the damping performance of an electrical power system against small disturbance(s), an 
excitation-based power system stabilizer (CPSS) is extensively used around the world. Although the CPSSs 
have been used widely for satisfactorily damping local-mode low-frequency oscillations, the outcome of 
CPSS may not be the best possible because of the intuitive nature of the tuning process and restrictive 
assumptions made. Later, the research was directed towards optimal control theory to develop an optimal 
state-variable feedback gain controller i.e. linear-quadratic-regultor (LQR) for stabilizing power systems 
against small perturbation(s). Consequently, reports [1]-[10] have appeared in the literature concerning the 
application of optimal LQR for stabilizing power systems. 

In the referred papers [1]-[10], the design of LQR is based on the following sequential process: 
i) The control (R) and state (Q) weights are chosen as diagonal matrices; ii) The state weighing matrix Q is 
assigned numerical values arbitrarily by an iterative procedure; iii) Optimal LQR is determined; and 
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iv) Closed-loop performance of an electrical power system equipped with optimal LQR is investigated. If the 
system performance is not satisfactory, the entire process has to be repeated until satisfactory damping 
performance is obtained. This trial-and-error (T&E) method doesn’t offer a systematic way of tuning Q; 
hence it is cumbersome, burdensome, and time-consuming. 

On the other hand, in recent decades, other researchers had been developed evolutionary algorithms 
to design optimal LQR [11]-[17]. In the referred papers [11]-[17], different evolutionary algorithms such as 
genetic algorithm (GA), particle swam optimization (PSO), big bang-big crunch (BB-BC), ant colony 
optimization (ACO), real-coded genetic algorithm, differential evolution (D.E.) algorithm, and Jaya 
algorithm are reported to design optimal LQR. And in the references [18]-[25], Gbest-guided artificial bee 
colony algorithm, backtracking search algorithm, adaptive backstepping approach, Whale optimization 
algorithm, an improved whale optimization technique, modified shuffled frog leaping (MSFL) algorithm, 
model reference self-tuning Takaji-Sugeno fractional-order proportional-integral-derivative (TSMFOPID) 
control technique and evolutionary programming based optimisation technique respectively for tuning CPSS. 
But, the evolutionary techniques are i) highly dependent on evolutionary algorithms to optimize and ii) 
computation time is more; hence their use is restricted for solving optimization problems. 

All these issues were addressed in this paper by proposing an analytical approach (AA) to tune 
LQR. The proposed technique of tuning LQR explores the correlation between the Lagrange-multiplier 
optimization technique and the algebraic-Riccati-equation (ARE). The proposed approach has the following 
advantages: i) It takes negligible time to tune LQR with the aid of derived algebraic equations; ii) It translates 
the performance objectives of the system from time-domain into cost-function; and iii) It enhances the 
robustness of the power system as Q varies in line with the operating condition. The proposed methodology 
explores its modular approach in tuning LQR. However, the design of a robust and decentralized LQR 
controller becomes impractical due to the non-availability of facilities to measure state variables (especially 
the rotor angle) in most multi-machine power systems (MMPS). This issue is addressed by considering the 
secondary voltage of transformer as a reference rather than infinite-bus voltage in modelling the power 
system [26]. 


2. MODELLING OF POWER SYSTEM 

An IEEE test system comprising 4-generators and 10-busbars shown in Figure | is chosen for the 
study. This section deals with the modelling and stability analysis of electrical power system. The modelling 
of power system is done based on modified version of Heffron-Phillips model [26]. This facilitates the 
proposed optimal LQR more practical. 


Figure 1. Single-line representation of 4-Generator and 10-Bus system 


The linearized dynamic equations of an i“’-generator of power system are; 


A6,, = (@p,ASm, — A8s,) (1) 

Mn = a (— G,,A5,, — Go,AEq, — Go, AVG, — DiASmn;) (2) 
i if AEq; 

AE;, = a (OEra = G4,Ads, == Gy, AVs; os aa (3) 
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; ee, r AEG; 
AEva, = fe (AVpss, — Gs, Abs, — Ge AEG, — Gv, AV) 7 a _ 


The state-space representation of (1)-(4) is; 


Xj = Ag x; + By ui + By, (5) 

Yi = Cy, Xj + Dg Ui (6) 
Since the coefficients G,, 1 Gy,» Gy, are negligibly small, the (5) and (6) are reduced to; 

Xj, = Ag,X; + Bg,u; (7) 

Vi = Cy, Xj + Dy Ui (8) 


The Laplace’s transformation of the state-space model of an i’-generator is; 


ye eee tev cet eeek, ot . 
sAEg,(s) = GaiTo, AE;,(s) To, Jost) + Tho, AE rq,(s) (9) 
sAé,,(s) = Aw;(s) (10) 
s4w;(s) = {Fr AT mn (S) = Fr [G1.485,(5) + G,,AEt(s)] - Fp DiAw(s)} (11) 
Kg. Kg. 
sAEpa,(s) = {FE ABCs) - 7g, LG5455,(8) + G¢,AEt,(s)] + AV er (8)} (12) 


Under the assumption that, AV,, f= 0; simplification of (9) and (12) yields (13); 


[G.,(1+sTE,) + Kz,G5,|%3; 


[(1+563,7do,)(1+5TE;) + KE; G6;6a,| 


AE;,(s) = A6d(s) (13) 


The Laplace’s transformation of electrical torque of an i’’-generator is; 


ATe,(S) = [G1,465,(s) + G2,4Ey,(s)] (14) 


By substituting (13) in (14) yields (15), 


[G4,(1 + s75;) + KE{Gs;|G,02, 


i 
1+ 563,740, )(2 + sTp,) + KE,G6,G3, 


AT. ;(S) ia {6 (( 28,0) = G(s) * A6,,(s) (15) 


3. PROBLEM STATEMENT AND ANALYTICAL APPROACH OF DESIGNING LOR 
This section explores the analytical approach of tuning state-weighing matrix (Q) associated with LQR. 
The feedback control law to minimize the cost function ‘J(u)’ of linear-quadratic-regulator is given by (16). 


u (t) = —Kx(t) (16) 
The cost function of LQR is given by (17). 

Ju) = fo x70). Q.x() + uO). R. u(t) dt (17) 
The Lagrange-multiplier optimization technique for optimizing linear state feedback gain vector (K) is, 

K=R 4B P (18) 
where, P is the matrix solution to the following reduced algebraic-Reccati-equation (ARE). 

A™P + PA+Q—PBR‘B'P =0 (19) 
The canonical form of electrical power system that is controllable [27] is given by; 
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x= Ax+Bu (20) 
y =Cx+Du (21) 
where, 
0 1 0 0 
we[Q ois 7 


B=[0 0 0 Balt (23) 
C=[0 Cy 0 Of (24) 
D = [0] (25) 


The (19) in accordance with controllable canonical form of power system becomes, 
ATP + PA+Q-—PBR-1B™P =0 (26) 


where Q, R and P are defined as (27)-(29). 


a 0 0 0 
_|9 gq 0 0 

C=l0 0 gs 0 (27) 
0 00 @ 

R=[r] (28) 


Pi1 P12 Piz) P14 


_ {P12 P22 P23 P2a 
Re Pi3 P23 P33 P34 ae 


P14 P24 P34 Pas 


The (18) in accordance with controllable canonical form of power system becomes, 

K =R-1BTp= “(Pus P24 P34 Daa | (30) 
For closed-loop systems, the actual characteristic equation is given by (31), 

|s1 —- A + BK|=0 (31) 


By substituting (22), (23) and (30) in (31) yields (32), 
2 2 2 2 
(s¢ + 59 (Chit — yy) +52 (Ha — yg) + s(t) + (EY - Aa) |= 0 02 


For 4"-order systems, the desired characteristic equation is (33), 
(s* + 4€@,s3 + 2w2(2&? + 1)s? + 4éw3s + wt) =0 (33) 


The 4"-row elements of matrix P are given by the following (34)-(37); 


Pia = (Wa + Aga) or (34) 
41 

Poa = (450 + Aga) oe (35) 

Daa = [2w2(2E? + 1) + Aas] 5 (36) 

Pag = (460, + Aaa) ar (37) 
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with the aid of ARE and (34)-(37) following equations are derived to get optimal Q. 


= gz (wh — 44) (38) 
Os py l2Aadas — Ada + 409(26? — 1] (39) 
B= oy l2ldn + Aarden + 3a% + Beng?(E? — 1)] — Ads} (40) 
Mt = oe [4a (28? — 1) - 2Ags — A44] (41) 


The (38)-(41) are the outcome of proposed analytical approach for designing an optimal LQR for specified 
natural frequency (w,), damping ratio(€), and scalar quantity R. 


4. DESCRIPTION OF POWER SYSTEM 

The loads of the [EEE test system considered are modelled as constant impedances. The total 
connected load in the system is 2734 MW. The machine data, line data, load flow data, and automatic voltage 
regulator (AVR) & excitation-system data of the test system are taken from [28]. The symmetric base system 
consists of two identical areas connected through a relatively weak tie line. Each area comprises two 
generating units with equal power outputs. The electro-mechanical modes of oscillation present in the system 
are modes present with in the plant (inter-plant modes) and low frequency inter-area mode. 


5. RESULTS AND DISCUSSION 

Here, a comprehensive discussion is made on results obtained against small disturbances i) step 
change in the voltage, ii) step change in the electrical torque under different loading conditions nominal load, 
light load, and heavy load. The comparative damping performances of generators that are fitted with either 
LQR tuned through analytical approach or LQR tuned via trial and error method or CPSS are shown in 
Figures 2-4. 


5.1. Case1: Nominal load 

In this case, all the generators are assumed to be loaded with their respective nominal loads. The 
Figure 2 shows the variation in slip (Sin) for 10% step increase in torque at generator Gy with tie-line power 
flow of 400 MW. It is observed that, the generator Gz equipped with proposed LQR exhibits much superior 
damping performance than its counter parts. The same is in the case for generators G;, G2, and G3. The 
settling times of generators noted from Figure 2 are tabulated in Table 1. 
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Figure 2. Variation in slip (S;,) when generator G4 operating at nominal load is subjected to 10% increase in 
torque with tie line power flow of 400 MW 


An analytical approach for LOR design for improving damping performance of ... (Sreenivas Uravakonda) 


56 Oo ISSN: 2502-4752 


Table 1. Settling time: nominal load: 10% step increase in torque at generator G4 
Settling time (Sec.) 


PaLDe CPSS _LQR tuned via T&E method — Proposed LQR tuned via AA 
G 8.6 8.6 2.6 
G 8.6 8.6 ~~ 
G 8.6 8.6 Ss 
Ge 8.6 8.6 zs 


5.2. Case2: Light load 

Here, all the generators are assumed to be loaded to 80% of their respective nominal loads. The 
Figure 3 shows the variation in slip (S,,) for 10% step decrease in reference voltage at generator G; with tie- 
line power flow of 100 MW. The generator G; equipped with proposed LQR settles at 1.9 Sec. with 
negligible overshoot and no undershoots after the disturbance. The same is in the case for generators G2, G3, 
and G4. The settling times of generators noted from Figure 3 are tabulated in Table 2. 


x 10° Generator-1 x10” Generator-2 
T T I I T T I I 
1 | ery Pdprctececcecee eee Pe eae CPSS a 1 | eee: ee oe CPSS a 
a : : eugines LQR-PSS tuned via T&E H He (pear LQR-PSS tuned via T&E 
= i LY fy. ix Proposed LQR-PSS a a. the : Proposed LQR-PSS 
2 0 aa —- — ne 2 sees : 
e€ uy Nol ' £ ' 
7) 7) 
i | eisascessaniactcdneontncisnetaevencisaaxe net cosuamaiasmenaxccesananees | ll iced aia panacea naam anaes auamesinaeh sRsmciaueeaeawaxe cnakes 4 
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Figure 3. Variation in slip (Sm) when generator G; operating at 80% of nominal load is subjected to 10% 
decrease in voltage with tie line power flow of 100 MW 


Table 2. Settling time: light load: 10% step decrease in reference voltage at generator G; 
Settling time (Sec.) 


pape CPSS _ LQR tuned via T&E method _ Proposed LQR tuned via AA 
G 4.2 10.4 19 
G 4.2 10.4 19 
G; 4.2 10.4 1.9 
Gs 4.2 10.4 19 


5.3. Case3: Heavy load 

In this case, all the generators are assumed to be loaded to 120% of their respective nominal loads. 
The Figure 4 shows the variation in slip (S,,) for 10% step increase in reference voltage at generator G3 with 
tie line power flow of 400 MW. In this case; after the disturbance, the generator G3 equipped with proposed 
LQR settles in 2.2 Sec. with no overshoot and negligible undershoot, whereas it is 10.6 Sec. and 8.8 Sec. 
when G3 is equipped with CPSS and LQR tuned via trial-and error method respectively. The settling times of 
generators noted from Figure 4 are tabulated in Table 3. 


Table 3. Settling time: heavy load: 10% step increase in reference voltage at generator G3 
Settling time (Sec.) 


Gen Ne: CPSS __LQR tuned via TRE method — Proposed LQR tuned via AA 
Gq 10.6 8.8 2.2 
Go 10.6 8.8 2.2 
G3 10.6 8.8 2.2 
Gy 10.6 8.8 2.2 
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Figure 4. Variation in slip (Sin) when generator G3 operating at 120% of nominal load is subjected to 10% 


6. 


increase in voltage with tie line power flow of 400 MW 


CONCLUSION 
In contrast with the LQR tuning via the trial-and-error approach, the LQR tuning via the proposed 


analytical approach leads to enhanced robustness of the power system as the state weighting matrix Q varies 
in line with the operating condition in the proposed methodology. The simulation results proved the 
superiority of the proposed optimal LQR in damping low-frequency electromechanical oscillations. It is 
observed that the proposed LQR improved the system dynamics by reducing the settling time and 
overshoots/undershoots. 
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